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Photonic qubits constitute a leading platform to disruptive quantum technologies due to their 
unique low-noise properties. The cost of the photonic approach is the non-deterministic nature of 
many of the processes, including single-photon generation, which arises from parametric sources 
and negligible interaction between photons. Active temporal multiplexing—repeating a generation 
process in time and rerouting to single modes using an optical switching network—is a promising 
approach to overcome this challenge and will likely be essential for large-scale applications with 
greatly reduced resource complexity and system sizes. Requirements include the precise synchro¬ 
nization of a system of low-loss switches, delay lines, fast photon detectors, and feed-forward. Here 
we demonstrate temporal multiplexing of 8 ‘bins’ from a double-passed heralded photon source and 
observe an increase in the heralding and heralded photon rates. This system points the way to har¬ 
nessing temporal multiplexing in quantum technologies, from single-photon sources to large-scale 
computation. 


I. INTRODUCTION 

Preparation and manipulation of exotic quantum 
states of light are at the heart of quantum informa¬ 
tion science and technology [T|. A central challenge for 
photonic systems is the non-deterministic nature of the 
generation of such states, which arises when consider¬ 
ing parametric photon sources [2], such as those based 
on spontaneous parametric down-conversion (SPDC) or 
four-wave mixing (FWM), and from the negligible inter¬ 
action between photons [3]- Parametric sources of single 
photons have been the workhorse for proofs of princi¬ 
ple to date because they generate photons in pairs, en¬ 
abling heralding in highly pure spatio-temporal-spectral 
modes [2]. However, parametric sources have a theo¬ 
retical maximum heralding efficiency of 25% [4], suffi¬ 
cient for some communication and sensing applications, 
but short of the best known threshold for computation 
(5] . As with other non-deterministic generation processes 
with heralded success signals, including fusion gates, for 
large-scale cluster states mm and all-photonic quantum 
communication [3, and ballistic entangled state genera¬ 
tion [MID], the success probabilities must be increased 
above relevant practical thresholds. 

A promising approach is to actively multiplex (MUX) 
these processes by operating several copies in parallel, 
such that the probability of at least one succeeding is 
high, followed by a low-loss switching network to route a 
successful output into the downstream system [B1IIIHI3]. 
Spatial multiplexing of heralded photon sources Eiini cur 
m, for example, has been successfully implemented with 
up to four heralded photon sources ushu]- Temporal 
multiplexing (19H2T] (see Fig. |T|) would enable repeated 
use of the same physical process, reducing resources, sys¬ 
tem size and indistinguishability requirements, at the 
cost of introducing delay lines and reducing the system 
clock rate. Temporal multiplexing has been proposed for 
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Figure 1: Illustration of the principle of temporal mul¬ 
tiplexing. A non-deterministic generation process is re¬ 
peated in time with period T; on heralded success, an active 
optical switching network and delay lines offset photons into 
output time bins spaced by an integer multiple of the input 
period and in sync with the system clock cycle. With a suffi¬ 
ciently low-loss switching network, the generation probability 
per clock cycle is increased. 


single-photon [IS EDI El [23] and entangled state gen¬ 
eration [S US El] [22, as well as for photon memories 
m, boson sampling schemes [25], and universal quan¬ 
tum computation [25], 

Here, we demonstrate active temporal multiplexing 
and use it to improve the success probability of a heralded 
single-photon source. By combining temporal with spa¬ 
tial multiplexing using a double-passed heralded source, 
only a single physical source was used to enable hybrid 
spatial-temporal multiplexing of eight effective source 
repetitions. We show active temporal multiplexing of a 
source of periodic photons (previous demonstrations have 
shown passive temporal multiplexing m and active stor¬ 
age, but not multiplexing, of photons from non-periodic 
sources using cavities [28] [29]) to increase the heralded 
photon rate, for a fixed noise level, by up to 76% com¬ 
pared to the same source without multiplexing. 
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Figure 2: Experimental set-up. Pulses from a femtosecond laser are upconverted using an LBO crystal, split into four copies 
using free-space delay lines, and passed twice through a PPLN crystal for down-conversion. Following separation of the photon 
pairs and pump using filtering, the heralding signals are analyzed by an oversampling FPGA while the signal photons are stored 
in long fiber delays. The FPGA configures the switching network to deliver the generated signal photons into a single spatial 
and temporal mode. PI and P2 indicate Pass 1 and Pass 2. 


Results 

Principle of operation 

In a heralded single-photon source (HSPS), a pulsed 
laser pumps a nonlinear material, spontaneously gener¬ 
ating photon pairs, called signal and idler photons, in a 
fixed time bin. Assuming spectral disentanglement m, 
the state after passing through the nonlinear material is 
given by an infinite superposition of Fock state pairs m 

IN = V 1 -!^ 2 (lN 1°), + E ? I> Ns) . (!) 

where i and s are the idler and signal modes and £ is the 
squeezing parameter determined by the pump power and 
the strength of the nonlinearity. Multi-photon pairs gen¬ 
erally result in detrimental effects in quantum circuits, 
necessitating low squeezing parameters so that the single¬ 
pair term in (1) dominates. 

Filters are used to separate the signal and idler photon 
and the pump, and a single-photon detector placed on the 
idler arm is used to herald the presence of the signal pho¬ 
ton. Under ideal conditions and with number-resolving 
detectors, the theoretical maximum single-photon emis¬ 
sion probability of a HSPS is limited to 25% [4], due 
to the presence of multi-photon pair terms in equation 
(1). While this single-photon emission probability is suf¬ 
ficient for small-scale quantum optics experiments, her¬ 
alded sources by themselves are not sufficient for scalable 
quantum technology m- 

A temporal multiplexing technique (also referred to as 
time multiplexing) [19, [20], which uses a HSPS, opti¬ 
cal switches, delay line loops, and electronics for feed¬ 
forward, can be used to boost the single-photon emission 
probability while keeping the multi-photon contamina¬ 
tion low (Fig. [l] shows the smallest example with one de¬ 


lay loop). In this scheme, the HSPS is pumped N times 
per clock cycle with laser pulses spaced by time T. Signal 
photons are stored in a long delay line buffer as detection 
signals from the idler arm are analyzed. When a single 
photon is heralded in one of the N time bins, a switch¬ 
ing network composed of delay line loops (with lengths 
of integer multiples of T) is driven into a configuration 
which offsets the photon into a single spatial-temporal 
mode. If multiple photons are heralded in several input 
bins, the switching network automatically discards the 
extra photons by moving them into adjacent bins, thus 
ensuring that only a single photon is output in the time 
bin in sync with the system clock cycle. 

With a sufficient number of time bins per clock cy¬ 
cle, a single-photon pair will be produced in at least one 
of the time bins with high probability. The probability 
of heralded single-photon emission from the multiplexed 
source is approximately (see Supplementary Material III 
for a detailed model): 

Vsingle ~ (l — 0- ~ Ptrig) ^ Psingle-, (2) 

where ptrig is the probability that the HSPS triggers 
during one time bin and p s ingie is the probability that 
the triggered emission is a single photon after passing 
through the lossy switching network M- With ideal op¬ 
eration and assuming a lossless switching network, 17 
heralded source repetitions enable a source with >99% 
single-photon emission probability [4], and assuming re¬ 
alistically small losses, ^8-16 heralded source repetitions 
enable a near-deterministic source with low multi-photon 
contamination for large-scale applications HU. Even 
when considering heralded sources operating with effi¬ 
ciencies far below the theoretical maximum, as is the case 
with all parametric sources demonstrated to date, mul¬ 
tiplexing can still be used to achieve an enhanced her- 
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Heralded Signal Photon (Coincidence) Rate 



Figure 3: Photon statistics from multiplexed and non- 
multiplexed sources, a, Triggering (idler singles) and b, 
and heralded signal photon (coincidence) rates vs. reference 
laser pump power for the 8x multiplexed, 4x multiplexed, 
and non-multiplexed sources. For clarity of presentation, 
data points from the non-multiplexed sources are shown as 
linearly interpolated region plots encompassing the range of 
data: blue (Pass 1, Delays 0-3) and orange (Pass 2, Delays 
0-3). "Delays 0-3" refers to the four effective non-multiplexed 
sources passively delayed in time. Theory lines for the multi¬ 
plexed sources are calculated from measured heralded source 
parameters, measured switch loss, and extrinsic loss effects. 

aided single-photon emission probability per clock cycle 
for a fixed multi-photon contamination probability, of¬ 
fering the possibility of new classes of experiments in the 
near-term. 

In this work we implemented a hybrid scheme using 
two spatially multiplexed sources fed into a temporal 
multiplexing set-up, doubling the number of effective 
source repetitions (Fig. |2j. By using a return pass in 
the opposite direction through the same nonlinear crys¬ 
tal, hybrid temporal and spatial multiplexing can be im¬ 
plemented with only a single parametric source. This 
scheme enables an additional enhancement to the single¬ 
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Figure 4: Heralded signal photon (coincidence) rate vs. 
CAR for multiplexed and non-multiplexed sources, a, 

Shows the full data set, and b, shows detail at low powers, 
where saturation effects due to electronics are small. "De¬ 
lays 0-3" refers to the four effective non-multiplexed sources 
passively delayed in time. Red points are for the 8x mul¬ 
tiplexed source, and blue points are for the non-multiplexed 
sources (Pass 1). Solid lines are based on a theory fit us¬ 
ing measured parameters. Dashed line shows a correction for 
extrinsic sources of loss based on the theory model. 


photon emission probability without an additional loss 
penalty on the generated photons, since the same depth 
of switches as the temporal multiplexing scheme is used. 

Implementation 

Our experimental set-up (Fig. [2]) uses a bulk 
periodically-poled lithium niobate (PPLN) down- 
conversion crystal phase-matched to produce idler 
photons at 671 nm and signal photons at 1547 nm from 
a pump laser at 468 nm. These wavelengths enable 
high-efficiency detection of the idler photon using silicon 
avalanche photodiodes (APDs) and low-loss transmission 
of the signal photon through switches and fiber delay 
lines. 
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Each pulse from a pump laser with a 80 MHz repetition 
rate (12.5 ns pulse spacing) is frequency doubled and split 
into four pulses spaced by ^3 ns using a series of free- 
space delay lines constructed from beam-splitters and 
mirrors. The four pulses then pass through the PPLN 
crystal and undergo collinear down-conversion, proba¬ 
bilistically creating photon pairs in four time bins in the 
first spatial mode (referred to as Pass 1); a return pass 
of the pump through the crystal, obtained by recycling 
the residual pump reflected from a mirror, creates four 
additional time bins in a separate spatial mode (Pass 2). 
For each pass, the four effective sources passively delayed 
in time are referred to as "Delays”, e.g. Pass 1, Delay 3. 
The spectra of the signal photons from different delays of 
the same pass were shown to have a high degree of sim¬ 
ilarity (mean of 97.9±1.8%, see Supplementary material 
II), and moderate similarity between the two different 
passes (mean of 92.8±5.2%, see Supplementary material 
II). Furthermore, an in-line polarizer was used to verify 
that the photons emitted from each source had identical 
polarization. 

The polarization-maintaining, active optical switching 
network (see Appendix A: Methods) is composed of low- 
loss fiber switches (^1 dB loss per switch, 500 kHz max¬ 
imum operation frequency) and two fiber delay loops 
matched to the free-space delay lines (see Fig. § m- 
Detection signal rising edges, which can fall in any of the 
four closely spaced time bins, are correctly discriminated 
using a fast oversampling Field Programmable Gate Ar¬ 
ray (FPGA) (see Appendix A: Methods), which then 
configures the switches for feed-forward multiplexing of 
the eight time and spatial modes. To avoid driving the 
switches faster than their maximum operation frequency, 
an asynchronous "idle time" of 2 fi s was programmed into 
the FPGA to limit the rate of detected heralding signals. 
During heralding detection, feed-forward processing, and 
switch configuration, the signal photons are stored in long 
delay lines of telecom fiber. Signal photons are detected 
using an InGaAs detector, gated from the idler photon 
detection events. 

Measured photon statistics from multiplexed and 
non-multiplexed sources 

Photon counting statistics were collected for the eight 
non-multiplexed sources and the multiplexed source. 
Triggering (idler singles) rates and heralded signal pho¬ 
ton (coincidence) rates (see Appendix A: Methods) are 
shown in Fig. [^(accidental rates are shown in Supplemen¬ 
tary material I). For all coincidence and accidental mea¬ 
surements, time-averaged dark counts, which had negli¬ 
gible effect at the measured rates, were subtracted from 
the totals. The data was taken for fixed "reference pow¬ 
ers", defined as the average power of the pump laser in 
front of the PPLN crystal used to generate the photons in 
the multiplexed source. Then, by blocking all beam paths 
in the free-space delay lines except for one at a time, a 
fraction of the reference power was used to pump each of 


the non-multiplexed sources individually (^25% for Pass 
1 delays and, due to power loss, ^12.5% for Pass 2 de¬ 
lays). This procedure allows for a fair comparison of the 
coincidence, accidental, and triggering rates between the 
multiplexed source and its constituent non-multiplexed 
sources. 

The triggering rates were affected by saturation loss 
caused by the large asynchronous "idle time" pro¬ 
grammed into the FPGA to avoid triggering detection 
faster than the limited switch operation frequency (500 
kHz). With the overall clock rate of the source set by the 
80 MHz repetition rate of the pump laser, the multiplexed 
source suffered from greater saturation effects at high 
powers due to this deadtime than the less-deterministic 
non-multiplexed sources, as can be seen in the sharp 
bending of the data away from a linear trend as the ref¬ 
erence power is increased in Fig. S*. 

The data was found to be in excellent agreement 
with our model of the non-multiplexed and multiplexed 
sources (see Supplementary material III), as shown in 
Fig. i According to the fit, the 4x time multiplexed 
source (composed of all Pass 1 delays) showed an in¬ 
crease of up to 175% in the triggering rate, for a fixed 
reference power, compared to the most efficient non- 
multiplexed source (Pass 1, Delay 3) and the 8x time 
and space multiplexed source (composed of both Passes, 
all delays) showed up to a 290% increase (see Fig. lo- 
At low reference powers (2.5-10.5 mW), the difference in 
coincidence rates for the 4x multiplexed source was not 
significant compared to the non-multiplexed source with 
the highest coincidence rates, while the 8x multiplexed 
source showed a higher coincidence rate than any indi¬ 
vidual non-multiplexed source, and a 17-69% higher rate 
compared to the mean from the non-multiplexed sources 
(Pass 1 only) (see Fig. [3 ]d). At high powers, the rates of 
coincidences for the multiplexed sources were suppressed 
due to the saturation of triggering events, as predicted 
by our model (Supplementary material III). 

Although the triggering and coincidence rate statis¬ 
tics provide evidence that our set-up implements tempo¬ 
ral multiplexing, the key measure of performance for our 
multiplexed source is the heralded photon rate for a fixed 
coincidence-to-accidental ratio (CAR). The CAR serves 
as a measure of noise due to single-photon emission and 
multi-photon contamination probabilities; this measure 
of noise cannot be inferred from coincidence rates alone. 
The heralded photon rates for fixed CAR for the 8x mul¬ 
tiplexed and non-multiplexed sources are plotted in Fig. 
[4] (data from the 4x multiplexed source is shown in Sup¬ 
plementary material I). For a fixed CAR, in the regime 
where saturation effects are small, the 4x multiplexed 
source did not show a significant increase in the heralded 
photon rate, and was limited mainly by the loss of the 
switching network (^4 dB loss for each path). However, 
the 8x multiplexed source exhibited an increased her¬ 
alded photon rate between 33-59%, for the same CAR, 
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over the best non-multiplexed source, and between 47- 
76% over the mean from the non-multiplexed sources 
(Pass 1 only), demonstrating a direct improvement (Fig. 
|3>). In our model we corrected for the effects of extrinsic 
loss on the rate of heralded photon production (dashed 
line Fig. [4^i, see Methods and Supplementary material 
III), indicating a potential improvement of ~114% for a 
wide range of CAR values compared to the expected her¬ 
alded photon rate from the non-multiplexed sources (also 
with extrinsic loss removed). This enhancement can be 
mapped to a comparable increase in the single-photon 
emission probability for a fixed multi-photon emission 
probability (see Supplementary material IV). 

Discussion 

We demonstrated temporal and spatial multiplexing 
of eight photon bins in a hybrid setup to enhance the 
heralded photon emission statistics compared to non- 
multiplexed sources. Although our demonstrated im¬ 
provement was limited by the maximum operation fre¬ 
quency of the switches, we note that even low rate, high 
efficiency multiplexed sources will likely find applications 
in the near term, due to the increased single-photon gen¬ 
eration probability per clock cycle. Therefore, a possible 
solution is to use a pulse picker to limit the repetition 
rate of the laser source to the maximum repetition rate 
of the switches. Ultimately, multiplexed sources with the 
highest single-photon emission rates will require the de¬ 
velopment of a lower-loss, high-speed optical switch (re¬ 
cent, promising prototypes include Kerr effect and 
electro-optic based |38j switches). 

Demonstrating Hong-Ou-Mandel interference between 
two independent multiplexed sources, which has not yet 
been demonstrated in any active multiplexing implemen¬ 
tation, is the next step in order to verify that photon pu¬ 
rity is preserved. In our multiplexed source, photons from 
different delays had similar spectral properties, identi¬ 
cal polarizations, and similar source efficiencies and cou¬ 
plings. Our setup of adjustable free-space delays allows 
for the fine tuning of the temporal delays of the photons 
to within a coherence length, indicating that photon in¬ 
terference between two of our multiplexed sources, up 
to the intrinsic limit of the PPLN sources themselves, 
should be possible. Improved multiplexing components 
and sources will further enable the temporal multiplex¬ 
ing of more complex generation processes, such as fusion 
gates mm or ballistic entangled state generation [SHE], 
where phase stability will be essential. 

Temporal multiplexing techniques will almost certainly 
be required in future large-scale quantum photonic cir¬ 
cuits in order to substantially reduce resource require¬ 
ments. Furthermore, hybrid temporal and spatial multi¬ 
plexing techniques will be important in order to optimize 
tradeoffs between spatial footprint and system clock rate. 


Integrated photonic components, including sources (e.g. 
[32]), switches (e.g. [S3]), filters (e.g. [34]), delay lines 
(e.g. m ), and detectors (e.g. [36]), are under develop¬ 
ment. Scaling down our set-up to a fully integrated pho¬ 
tonic chip with low-loss components will enable a tempo¬ 
ral multiplexing template capable of realizing new classes 
of quantum information experiments and technology. 


APPENDIX A: METHODS 
Experimental set-up 

A mode-locked, Ti:Sapphire laser ( M Tsunami M , Spectra 
Physics) produced ^150 femtosecond pulses at 936 nm; 
a LBO crystal (Newlight) was used to frequency convert 
to 468 nm. To enable low-loss, near 50-50 splitting of 
the pulsed pump beam, laser line non-polarizing beam¬ 
splitters (Newport) were used in the free space delays. 
The PPLN crystal (Covesion) was 3 mm long and phase- 
matched at 110° C using an oven and temperature con¬ 
troller. Dichroic mirrors (Semrock) were used to separate 
the signal and idler photons from the pump, and Pellin- 
Broca prisms were used for further spatial filtering. A 
bandpass filter centered at 671 nm (Semrock) was used 
on the idler arms of each pass for further filtering. 

Polarization-maintaining switches (Agiltron, ^1 dB 
loss per switch, 500 kHz max operation frequency) 
were based on an electro-optic material. Standard tele¬ 
com fiber was used for the long delay buffer (^200 m, 
^90% transmission) and polarization-maintaining, low- 
dispersion fiber (Corning) was used for the variable de¬ 
lay line loops (^0.65 m and 1.30 m, ^95% transmission). 
Fiber polarization controllers (Fiberpro) were used be¬ 
fore the polarization-maintaining switching network to 
match the polarizations of photons from the two passes. 
To enable reliable comparison between multiplexed and 
non-multiplexed source measurements, a MEMS switch 
with nearly balanced loss was used to route the photons 
into or around the multiplexing switch network. 

Photon detection 

Idler (triggering) photons were detected using silicon 
avalanche photodiodes (APDs) (PerkinElmer). Pump 
leakage and dark counts were found to be negligible on 
the idler arms. 

Idler photon detection signals were discriminated with 
an oversampling FPGA (Xilinx Spartan 6) using internal 
delays and a 80 MHz, M locked-to-clock M reference input 
from the Ti:Sapphire laser. The FPGA was designed so 
that after signal detection, an M idle time” of 2 /is became 
active to avoid further detection. For every detected sig¬ 
nal (regardless of time bin), a gating signal was output 
with a constant delay with reference to the original input 
clock, to ensure correct heralding of the temporally mul¬ 
tiplexed photons. The total (unoptimized) internal delay 
of the FPGA logic was ^60 ns. 

Signal photons were detected using InGaAs detectors 
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(ID Quantique). Coincidence counts, joint detection be¬ 
tween idler and signal photons from paired generation 
events, were collected using gated detection of the signal 
photon from idler detection signals from the FPGA. Acci¬ 
dental counts, joint detection between idler and gated sig¬ 
nal photons from unpaired generation events, were then 
collected by shifting the temporal delay of the FPGA 
input clock by a multiple of the clock cycle. Pump leak¬ 
age in the signal arms was found to be negligible at the 
measured powers. Dark counts detected by the InGaAs 
detectors in gated mode were measured by blocking the 
signal arm path; these time-averaged counts were then 
subtracted from the count totals. The gate width used 
was 1.8 ns. 

Extrinsic sources of loss 

Extrinsic source of loss in the setup include: 1) Loss 
due to measurement apparatus. A small amount of extra 
loss (4%) on the multiplexed source was due to asym¬ 
metric loss of the MEMS switch used to switch between 
multiplexing and non-multiplexing channels for measure¬ 
ment. 2) Loss due to the deadtime of available electronic 
amplifiers. The two electronic amplifiers used to amplify 
the signal from the APD and into the FPGA have dead- 
times of ~ 0.1 ii s, resulting in missed pulses from the 
APD. In principle, much faster amplifiers with negligible 
deadtimes can be used to eliminate this source of loss. 3) 
Loss due to the limited switch repetition rate (500 kHz). 
An asynchronous M idle time" of 2 /is was programmed 
into the FPGA to avoid driving the switches faster than 
their 500 kHz maximum operation frequency. The switch 
repetition rate is set by the switch driver board; the 
switches themselves have a faster intrinsic rise and fall 
time of 300 ns (~ 3 MHz). 
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Supplementary Material: Active Temporal Multiplexing of Photons 

II. ACCIDENTALS AND CAR PLOTS 

A plot of accidental rates against reference power for the multiplexed and non-multiplexed sources is shown in Fig. 
Sla. A plot of coincidence rates against CAR, including data from the 4x multiplexed source and from Pass 2 is 
shown in Fig. Sib. 


III. SOURCE DETAILS AND SPECTRAL DATA 


Although Pass 2 should ideally be identical in operation to Pass 1, it differs from Pass 1 in two aspects. First, 
the laser pump power passing through the down-conversion crystal has been reduced by ^50% compared to the first 
pass, due to loss from the crystal without anti-reflect ion coating (^15% pump loss per facet) and additional loss from 
dichroic mirrors used for filtering. Second, Pass 2 is affected by unwanted “back-reflected” photons from the Pass 
1, again due to a lack of anti-reflect ion coating and the cut of the crystal. These back-reflected photons result in 
false heralding events, degrading the quality of the multiplexed source. By optimizing the position of the crystal, 
the back-reflected photon contamination in Pass 2 can be minimized (^20% of the total counts on the idler arm and 
^ 8 % of the total counts on signal arm). This results in a slight change in the CAR, although the larger effect is 
additional saturation loss (described in the main text) due to the increased triggering rate. These features were taken 
into account during the modelling and analysis of results (Supplementary material IV). 

Spectral data from all eight effective sources was collected by detecting signal photons using InGaAs detectors after 
passing through a tunable filter (Yenista), with time-averaged dark counts subtracted (Fig. S 2 ). To quantify the 
similarity between photons from different passes and delays, upper bounds on the indistinguishability parameters 7 
|S39| between each pair of sources were calculated using integral overlaps of the fit spectral amplitude functions (Table 
[I]). The integral overlap formula is 




Figure SI: Accidental and CAR plots, a, Accidental rates vs. reference laser pump power for the 8x multiplexed (red), 
4x multiplexed (green), and non-multiplexed sources. For clarity, data points from the non-multiplexed sources are shown as 
linearly interpolated region plots encompassing the range of data: blue (Pass 1, Delays 0-3) and orange (Pass 2, Delays 0-3). 
PI DO indicates Pass 1, Delay 0, and similarly for the other labels. Theory lines for the multiplexed sources are calculated from 
measured heralded source parameters, measured switch loss, and extrinsic loss effects. Inset shows detail at low powers, where 
saturation effects due to electronics are small, b, Heralded signal photon (coincidence) rate vs. CAR for multiplexed 
and non-multiplexed sources. Red points are for the 8x multiplexed source, green points are for the 4x multiplexed source, 
blue points are for the non-multiplexed sources (Pass 1), and orange points are for the non-multiplexed sources (Pass 2). PI DO 
indicates Pass 1, Delay 0, and similarly for the other labels. Solid lines are based on a theory fit using measured parameters. 
Dashed line shows a correction for extrinsic sources of loss based on the theory model. 
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Figure S2: Source spectral characterization. Normalized spectral data from signal photons from both passes. DO indicates 
Delay 0, and similarly for the other labels. A max is mean peak position and FWHM is mean full-width half maximum for all 
delays in a pass. Curves are best fits to a Gaussian. 


Pass 1 

Delay 0 

Delay 1 

Delay 2 

Delay 3 

Delay 0 

1.0 

0.966 

0.988 

0.952 

Delay 1 


1.0 

0.994 

0.997 

Delay 2 



1.0 

0.986 

Delay 3 




1.0 


Pass 2 

Delay 0 

Delay 1 

Delay 2 

Delay 3 

Delay 0 

1.0 

0.991 

0.946 

0.997 

Delay 1 


1.0 

0.979 

0.993 

Delay 2 



1.0 

0.960 

Delay 3 




1.0 



PI DO 

PI D1 

PI D2 

PI D3 

P2 DO 

0.997 

0.979 

0.925 

0.993 

P2 D1 

0.947 

0.903 

0.828 

0.941 

P2 D2 

0.967 

0.936 

0.877 

0.969 

P2 D3 

0.938 

0.893 

0.819 

0.934 


Table I: Tables showing upper bounds on indistinguishability parameters between photons from different passes 
and delays. Top table shows upper bounds on indistinguishability parameters from photons from different delays in Pass 1. 
Middle table shows upper bounds on indistinguishability parameters from photons from different delays in Pass 2. Bottom 
table shows upper bounds on indistinguishability parameters from photons from delays from different passes. "P" and "D" 
refer to the pass and delay number, respectively. 







































3 


7 = 



A M ^B M dw 


2 


(SI) 


where ipA (w) and 'ipB (<^) are the spectral wavefunctions of the photons from two different effective sources. 7 obeys 
0 < 7 < 1, where 7 = 1 corresponds to perfect indistinguishability and 7 = 0 corresponds to perfect distinguishability. 
We emphasize that this is only an upper bound on indistinguishability, and does not consider other degrees of 
freedom or joint spectral entanglement between signal and idler photons. Between different delays from Pass 1, the 
mean overlap was 0.980T.017, and between different delays from Pass 2, the mean overlap was 0.978T.020. This 
similarity can likely be further improved by passing the pump beam through a single mode filter such as a pin hole 
setup before the PPLN crystal, eliminating possible angular and intensity deviations in the pump beams for different 
delays. Between delays from different passes, the mean overlap was 0.928T.052. The difference in spectra can likely be 
explained by slight deviations in the Pellin-Broca filtering and coupling setups for the different passes, a feature that 
can affect any two independent sources and is not inherent to the multiplexing setup. Careful alignment optimizations, 
or further narrowband filtering, can help eliminate this source of distinguishability. 


IV. THEORETICAL MODEL 
A. Heralded single-photon source (HSPS) Model 

We assume that strong spectral filtering on the idler arm and moderate filtering on the signal arm are sufficient 
such that the states produced by our sources are close to an idealized two-mode squeezed state with disentangled joint 
spectra. Each source is characterized by three parameters: 1) an input power seed P see d which relates the amount of 
input pump power to a reference photon generation probability from the crystal (p pa i r = 0.1). 2) Idler transmission 
rji, which includes all sources of loss on the idler photon, such as from filters, coupling, and detector inefficiencies. The 
effect of nonlinear loss due to saturation effects from electronics and deadtime are treated in the next subsections. 3) 
Signal transmission r] s which includes all sources of loss on the signal photon. 

The relation between pump input power and effective squeezing parameter £ for a down-conversion source is given 

by 


£ = tanh ^ c\[P ^ , (S2) 

where P is the input power in units of mW and c is a coupling constant in units of mW -1 / 2 . The squeezing 
parameter £ see d corresponding to p pa ir = 0.1 is £ 8 eed ~ .335715, allowing for the extraction of the coupling constant 
c given a seed power P see d • With this coupling constant, the squeezing parameters corresponding to any input power 
can be found. 

Given a squeezing parameter £, idler transmission rfr , and signal transmission rj s , the probability for the idler 
(herald) arm to trigger using threshold detectors is given by [S14 ]: 


\£vi 

1 -KI 2 (!-%)■ 


(S3) 


Given that the herald has triggered, the probability for a single-photon emission to lead to a detection event on the 
signal arm is 


Psingle 


(i-ICI 2 ) 


Vs 


1 - 


(l£| 2 (l -Vs)) (1 -Vi) 


1 - I€I 2 (1 -Vi) 


1 - i^i 2 (i - ^)i 2 b - ier (i - (i - ^)l 2 


(S4) 


Given that the herald has triggered, the probability for a multi-photon emission to lead to a detection event on the 
signal arm is 
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Z T d = (1 — 




P multi 


Z T D 


Pi 


Psingle 


rig 


(! - Vi) 


1 , (1 - rjs) (1 ~ Vi) _ 

1-|^| 2 1-ICI 2 (1- Vs )(l- Vi ) 1- ICI 2 (1-77,) 


0 --Vs) 


i- l£l 


(S5) 


(S6) 


Using these expressions we can find the probability of coincidence and accidental detection, and CAR (coincidence- 
to-accidental ratio). 


B. Coincidences, accidentals, and CAR 

We can find the probability of coincidence detection using threshold detectors, assuming the contribution from dark 
counts has been subtracted and pump leakage is negligible. 

The probability of a coincidence is given by 


Pc Ptrig (Psingle "b Pmulti) ‘ 


When seeded with a pump laser with repetition rate R, the expected coincidence rate is 


(S7) 


R c = Rx Pc- (S8) 

An accidental occurs when a coincidence between an idler and signal photon generated from two different pump 
pulses occurs, assuming the subtraction of dark counts and negligible pump leakage. 

Without triggering from a herald, the signal photon trigger probability is given by 


Ptrig 


I£|V 


i-ira-%) 

Then a good approximation to the accidental probability is given by 


(S9) 


P A = PlrigPtrig- (S10) 

The rate of accidental detection events is then 

Ra = Rx Pa- (Sll) 

The expected CAR, coincidence to accidental ratio, is then 

CAR = R c /R a . (S12) 


C. Electronics saturation 

In practice, saturation effects due to electronics and detector deadtimes will further affect the source performance. 
Assuming a detector has a deadtime of d, and that events are approximately uniformly distributed, the true rate of 
counts T from a detected rate of counts D can be approximated as [340] 

D 

1 -Dd' 


T = 


(S13) 
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Source 

Vi 

Vs 

P see d (mW) 

R 2 

Pass 1 Delay 0 

0.015 

0.0019 

5.2 

0.993 

Pass 1 Delay 1 

0.015 

0.0019 

6.8 

0.990 

Pass 1 Delay 2 

0.016 

0.0021 

5.6 

0.988 

Pass 1 Delay 3 

0.017 

0.0018 

4.6 

0.993 


Table II: Table showing source parameters for Pass 1 . Source parameters were determined using a numerical optimization 
fit to the model. R 2 is the mean R 2 from the triggering rate, coincidence, and accidental fits. 


Conversely, given an expected rate of counts T, the detected rate of counts will be 


T 

dT + l’ 


(S14) 


By applying these equations several times, the effect of several electronic deadtimes in series can be modelled. Most 
of the saturation effects in the experiment affect the idler arm, so these equations can be used to find an effective 
Ptrig probability. Then, an effective ry which includes loss due to saturation can be determined using equation 


S3 


D. HSPS Model fitting to data: Pass 1 

Pass 1 has four time bins passing through the down-conversion crystal, which we will call Delay 0, Delay 1, Delay 2, 
and Delay 3. We will label these sources with the tuple ( pass , delay ) corresponding to the pass and delay associated 
with the state source, for example Source (1,3) identifies the source from Pass 1 Delay 3. 

Although in the ideal case each source would be exactly the same, in practice each source is slightly different: due 
to slight imperfections in alignment and optics, each source may feature different relative pump input powers, pump 
coupling, and different signal and idler loss rates. 

Using the above model of the HSPS, a numerical optimization was performed to find the optimal input power seed 
Pseed') idler transmission value ry, and signal transmission value rj s to match recorded values of heralding triggers, co¬ 
incidences, and accidentals. We will label these parameters for each source with the corresponding tuple: ^P ass ^deiay) , 
^(pass,delay) ^ p(pass,delay). moc [ e i s w hich maximized the mean R 2 statistic for trigger counts, coincidences, and 

accidentals were found for all delays. These are shown in Table [TT| 


E. Time Multiplexed Source Model: Pass 1 

The time multiplexed source is pumped by laser pulses which have been split into four pulses using free-space delay 
lines (Main text: Fig. 2), so the time multiplexed source model is constructed from the four HSPSs models from 
section D. Each HSPS in the time multiplexed source is pumped with a fraction of the total power used to pump the 
multiplexed source; the experimentally measured values were 0.2375, 0.2693, 0.2258, 0.2586 for Delay 0,1,2, and 3, 
respectively. The model for the time multiplexed source follows that of [SI 4 ], slightly modified to account for each 
time bin source being slightly different. In this section we neglect saturation effects; the same methods from section 
C can then be applied to the expected detection events to derive detection events with saturation. 

The probability for the multiplexed source to trigger is given by 


MUX,1 

Ptrig 


= 1 - 


( X -Plrig 0) ) (l ~ Ptrig^) ~ Ptrif) • 


When seeded with a pump laser with repetition rate R, the heralding rate is then 


(S15) 


r>MUX,l 
^trig 




(S16) 


Let ?^ as5 ’ deZay) correspond to the transmission due to the switching network and variable delay lines for source 
( pass , delay). Then the probability the heralded state is a single photon emission from source ( pass , delay), including 
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the effect of the switching network with variable delay lines is p S s ingie S ’ delay \ given by replacing each instance of r] s 
in equation S4 with r] s r]i^ ss,delay \ Similarly, the probability the heralded state is multi-photon emission from source 
(pass, delay) is p^m , is given by replacing each instance of rj s m equation S5 with p s Psw • 

The coincidence probability is 


jjMUX, 1 _ p' 


(bo)' 


+ (i-iST) (4 Mr + (i -riS 1 ’) (rg* + (i-^T)if w )). 


(SIT) 


where P, 


(1 ,delay)' i,(l,delay) 

~ Ptr 


s (1 ,delay) . s (1 ,delay) \ 
P single ' Pmulti 1 


C "trig single ' "multi J * 

When seeded with a pump laser with repetition rate R, the expected coincidence rate is 


r>MUX, 1 

K c 


= R x 


tjMUX,1 


(S18) 


The accidental probability is 


tjMUX 

VA 


,1 


z,(1,0) s',(1,0) 
"trig "trig ' 


MS 0) )( 


PtHg^Ptrig’^ + 


V’ (M) 

ring 


){A, 


(1,2) s', (1,2) 
rig "trig 


The rate of accidental detection events is then 


+ 


MS 21 ) 


i,(W) 

ring 


^'.(1.3) 

ring 


))• 


(S19) 


R MUX,l =Rx pMUX, 1 _ (S20) 

The expected CAR, coincidence-to-accidental ratio, is then 

CAR mux + = R^ ux /R^ ux . (S21) 

Using measured values of switching loss and electronics saturation, the model was found to be in good agreement 
with the data (Main text: Fig. 3). 


F. HSPS Model fitting to data: Pass 2 


Pass 2 has two features different from Pass 1 (Supplementary material III). We used a simple model describing 
the effect of the back-reflected photon contamination. We assume the back-reflected photons on the signal arm are 
negligible. On the idler arm, we assume that every delay has the same probability of triggering due to a back-reflected 
photon as a fraction / of the “true” triggering probability. 


i,(2,delay) _ i,(2,delay) n 

Ptrigback Ptrigtrue J‘ 


(S22) 


Then the probability of “correctly” triggering is 


i,(2,delay) _ i,(2,delay )_ i, (2, delay) x i,(2,delay) ,i, (2, delay) x 

Ptrig, correct Ptrig true V 1 Ptrigback / ' Ptrig true \Ptrigback /’ 

and the probability of “incorrectly” triggering is 


(S23) 


i, (2,delay) _ i, (2, delay) (_ i,(2,delay) \ 

Ptrig,incorrect Ptrigback l Ptrigtrue J ’ 

Then the total probability for one arm to trigger is 

i,(2,delay) _ i,(2,delay) i,(2,delay) 

Ptrig Ptrig,correct ' Ptrig,incorrect' 


(S24) 


When seeded with a pump laser with repetition rate R, the heralding rate is then 
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Source 

Vi 

Vs 

P see d (mW) 

R 2 

Pass 2 Delay 0 

0.018 

0.0024 

6.3 

0.971 

Pass 2 Delay 1 

0.017 

0.0021 

6.7 

0.982 

Pass 2 Delay 2 

0.016 

0.0023 

6.8 

0.978 

Pass 2 Delay 3 

0.015 

0.0020 

6.9 

0.987 


Table III: Table showing source parameters for Pass 2. Source parameters were determined using a numerical optimization 
fit to the model. R 2 is the mean R 2 from the triggering rate, coincidence, and accidental fits. 


7 ~>i,(2 ,delay) 
U trig 


R x p\ 


, (2, delay) 
rig 


(S25) 


Following the techniques in m, we can derive the probability of single and multi photon detection given that the 
heralding detector has not triggered due to a paired idler photon: 


i, (2, delay) 

Psingle no trig 


(i- l€l 2 (i — ^?i)) Vs(i~Vi)\Z \ 2 


(S26) 


i, (2, delay) 

Pmulti no trig 


( i -|£| 2 ( i -%)) l £| 2 


(! ~ Vi) 

1 - |^| 2 ( 1 - 77 ,) 


(1 - rj a ) (1 - rjj ) \ 

1 -| tf(l-ris)(l-rk)) 


i, (2, delay) 

Psingle no trig' 


Then the coincidence probability is 


p(2, delay) 


_ b(2,< delay) ( i, (2, delay) i,(2,delay)\ 

" Ptrig,correct \P single ' Pmulti J 


i, (2, delay) 

U trig,incorrect 


i, (2, delay) 


i, (2, delay) 


Psingle no trig ' Pmulti no trig 


.) 


The expected coincidence rate is 


(S27) 


(S28) 


p(2 ,delay) 

K c 


= R x 


p(2, delay) 

r c 


(S29) 


The accidental probability is 


p(2, delay) i, (2, delay) s, (2, delay) 

U4 Ptrig Ptrig 


(S30) 


and the expected accidental rate is 


^>(2 ,delay) = ft x p(2, delay) 

The model which maximized the mean R 2 statistic for trigger counts, coincidences, and accidentals was found for 
all delays. These are shown in Table IH 

G. Time Multiplexed Source Model: Pass 2 


The probability for the MUX source from Pass 2 to trigger is 


MUX, 2 
Ptrig 


= 1 - (! -Plri 9 0) )(l -PtrT ) )( 1 -Ptrig 2) )(l -Ptrif), 


and the expected heralding rate is 


(S31) 
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^trig 
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R x p trig . 
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The coincidence probability is 
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where p^ delay ^> j s given by p& delay ) except with every instance of p s replaced with r] s r]iw delay \ 
The expected coincidence rate is 
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(S34) 


The accidental probability is 


pM[/X,2 
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rtrig "trig ~ 
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where p trig is given by p trig except with every instance o± r] s replaced with p s rjsw 

The expected accidental rate is 


i, (2,0) s',(2,1) 
Jrtrig Jrtrig 


))• 


(S35) 


rjMUX ,2 t-> w tjMUX ,2 

jA — it X -J /^ 


(S36) 


H. Time and Space Multiplexed Source Model: Both Passes 

The probability for the complete multiplexed source to trigger is 




and the expected heralding rate is 


The coincidence probability is 


R^ q x = Rxpf 


rig * 


pMUX _ pMUX 

c c 


,1 . MUX,l\ 

+ { l ~ Ptrig ) 


d MUX,2 
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The expected coincidence rate is 


The accidental probability is 


R% ux = Rx P™ ux . 


U MUX _ t)MUX, 1 . MUX,l\ U MUX,2 

' .1 “ 1 A I ^ 1 Ptrig ) R 4 


The expected accidental rate is 


RMUX = r X pMUX' 


The expected CAR is 


(537) 

(538) 

(539) 

(540) 

(541) 

(542) 


car mux 


pMUX i pMUX 


(S43) 


Using measured values of switching loss and electronics saturation, the model was found to be in good agreement 
with the data (Main Text: Fig. 4). 
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Figure S3: Extracted theoretical single-photon emission probability for multiplexed and non-multiplexed 
sources, a, Model with extrinsic sources of loss removed, b, Model with all sources of loss except switch loss and power loss 
removed. 


V. SINGLE-PHOTON EMISSION ENHANCEMENT 

Using the model with extrinsic sources of loss removed, we can extract the heralded single-photon emission prob¬ 
abilities for fixed multi-photon emission probabilities of the 8x and 4x multiplexed sources and the best performing 
non-multiplexed source (Fig. S3a). We include the loss of final measurement detector. 

To further examine the long-term prospects of the multiplexed source, we can also remove all sources of loss, except 
those due to the multiplexing switches and power loss affecting the second pass. These sources of loss include filtering, 
coupling, and detector inefficiencies (although we still assume the use of threshold detectors). The extracted single¬ 
photon emission probabilities for fixed multi-photon emission probabilities of the 8 x and 4 x multiplexed sources and 
the best non-multiplexed source are shown in Fig. S3b. 















